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Exact ground states of interacting electrons on the diamond Hubbard chain in a magnetic field 
are constructed which exhibit a wide range of properties such as flat-band ferromagnetism and 
correlation induced metallic, half-metallic or insulating behavior. The properties of these ground 
states can be tuned by changing the magnetic flux, local potentials, or electron density. 
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Condensed matter systems with macroscopic degen- 
eracies react very sensitively on internal or external per- 
turbations and thus give rise to fascinating emergent be- 
havior. Well-known examples are electrons in a magnetic 
field in two dimensions |l[ and spins on lattices with ge- 
ometric frustration Q|. Dispersionless ("flat") electronic 
bands in solids also lead to macroscopic degeneracies. Re- 
cent advances in nanotechnology allow for the possibil- 
ity to design simple structures, which have flat single- 
electron bands,i.e., when electronic correlations are ne- 
glected 0, Him, Hi. The controlled setup of optical lattices 
for cold atoms using standing wave laser light-fields also 
allows to realize a variety of lattice models of interacting 
fermions and bosons @, @] • The ability to generate flat 
bands by changing external parameters such as a gate 
voltage or a magnetic field, would make even a direct 
manipulation of macroscopically degenerate systems pos- 
sible. For example, the understanding of ferromagnetism 
in organic compounds [5[ could be improved, and the 
tuning of instabilities or the switching between different 
phases would permit direct applications in spintronics. 

Most exact flat-band results concern flat lowest-energy 
bands, and provide solutions for ground states or the 
low temperature thermodynamics [H, [t| . For such a case 
Mielke and Tasaki proved that ferromagnetism is stabi- 
lized at or near half- filling [loj] . Lieb's ferrimagnetism 
emerges on a bipartite lattice with a macroscopically de- 
generate energy level exactly in the middle of the spec- 
trum [11]. For the more general case, when the inter- 
action leads to an effective flat band above a dispersive 
band, exact results were derived for the periodic Ander- 
son model [13] , by exact diagonalization of small clusters 
and ah initio band-structure calculations ||, or analyti- 
cally for the two-particle problem [f| . 

In this Letter we concentrate on one of the simplest 
lattice electron models, the Hubbard model on a dia- 
mond chain, as found, e.g., in azurite. In the diamond 
chain flat bands occur already in the one-electron picture. 
External potentials or a magnetic field can give rise to 
additional flat bands. For selected parameter sets, and 
by appropriately tuning local potentials or the magnetic 
flux, we construct exact many-body ground states on the 



diamond chain, which are either insulating or conduct- 
ing and fully or partially spin polarized. Our results thus 
open new routes for the design of spin valve devices and 
gate induced ferromagnetism. 

Fig. 1 shows the diamond chain, whose sites are de- 
noted by i + r s , where i and r s (with sublattice index 
s = 1,2,3 and r3 = 0) denote the unit cells and the three 
sites inside, respectively. The Bravais vector of the chain 
is a; periodic boundary conditions along the chain are as- 
sumed. Given N electrons their density is n — N/(3N C ), 
where N c is the number of unit cells. The Hamiltonian 
for the diamond Hubbard chain has the form 
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where the kinetic and the interaction part are written 
in k space and position space, respectively. Here and in 
the following sums or products over k and i extend over 
the N c vectors enumerating the unit cells, cj a creates an 
electron with spin a at site j, the local density is given by 
riy a = c| )0 .Cj !(T , and c Si k,<r is the Fourier transformed sub- 
lattice operator. The elements of the symmetric matrix 
M SjS /(k) in the kinetic energy are 



M u = M 2 , 2 = e, M 3 , 3 = 2*|| cosafc, M x , 2 
M is=1 , 2 ), 3 = 2tcos[(ak + (-l) s 5)/2}. 
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FIG. 1: Diamond Hubbard chain; the cross-shaped region 
depicts a localized (Wannier) eigenfunction for t± = tii = 
and flux S = tt/2 (see text). 
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FIG. 2: Single-electron band structure: (a) t± = tn = 0, 



S = tt/2. (b) e = -tx+tZ A 5 : 
degenerate for txty = 1/4. (c) ty 



7r; the upper flat bands are 
= 0, et± cos 8 = t j 



Here, t, in denote the hopping matrix elements to 
nearest neighbors (NN) and to next NN sites (perpen- 
dicular and parallel to a), respectively. A Peierls phase 
factor exp(i(5/2) with 5 = 27r<i>/$o accounts for a perpen- 
dicular magnetic field B; $ is the flux threading the unit 
cell and $o = hc/e is the flux quantum. Choosing the 
vector potential A || a in the field dependent hopping 

amplitudes iy(B) = tj j/ (0) exp[(i27r/$ ) f- A • dl] the 
magnetic field does not alter t± and ty . The on-site po- 
tential e acts on the sublattices s = 1, 2. U is the on-site 
Hubbard repulsion. In the following energies will be in 
units of 2t. 

The diagonalization of the kinetic energy part in Eq. 
((T|) leads to a cubic eigenvalue equation for the band dis- 
persion E(k). With Pi = (1 + t\) + cos (5 cos ak and 
Po = t± [ cos <5 + cos afc + 1± (e — 2tii cos ak)] the dispersion 
J5 w (k) (i/ = 1, 2, 3) is determined by 

(E - e) 3 + {E-e) 2 (e- 2t\\ cos ak)-(E~e)P 1 = P . (3) 

Depending on the choice of tj_, i|| , the magnetic field, and 
the potential e, the single-electron band structure con- 
tains one, two or even three exactly flat bands (Fig. 2). 
It should be stressed that these bands lose their meaning 
as soon as the interaction is switched on. 

Vidal et al. recently presented a detailed study of two 
electrons on the diamond Hubbard chain in the limit 
tx = ill = e = Q. They showed that for half a flux 
quantum per unit square (6 — tt/2) the excited singlet 
eigenstates are localized if U = 0, but become delocal- 
ized for U > 0. Apparently, the interaction U is able to 
induce subtle correlations leading to conducting states, 
which led them Q to speculate that such a derealiza- 
tion also holds for a finite electron density. 

Here we construct exact many-particle ground states 
for quite general cases and thereby also resolve some of 
the issues raised in Q. The strategy for deducing exact 
ground states was described before in the context of the 
periodic Anderson model [l2j ]. but has not yet been ap- 
plied to Hubbard models. The key steps are to first cast 
the Hamiltonian in positive semidefinite form and then 
to construct an explicit eigenstate with minimal energy. 

Solution I. Flat-band ferromagnetism: We start with 



and S = tt/2 (" Aharonov-Bohm cage" Q), in which case 
Eq. ([3]) provides three flat single-electron bands with 
energies E 2 = e, E 2 ±i = (e =F \/e 2 + 4)/2 (see Fig. 2a). 
Introducing the canonical fcrmionic operators C v \^ 

1 6 (_) 
C2±i,i,a = -^[FxQ^ T 2F T c^] , C 2Xa = (4) 



with Fl = 1 T Qg° = Q£(T5) ± 0£ J (±£) 

and Qf =1 ^ 2) (^) = e- i5 / 2 c i+ri , CT + e+^c;^^, the 
Hamiltonian can be written in the form 

3 

fl = EE^A+^ ( 5 ) 



i(2) 



Since C\, i a C vX „ and Hjj in Eq. ([5]) are positive semidef- 
inite operators, the ground state of N < N c electrons at 
U > is obtained by filling up the eigenstates as 

N 

\* I g m = i[cl itai m (6) 

i=l 

|0) denotes the vacuum. The ground-state energy is 
Eg = EsN. The operators & u ; a create electrons in local- 
ized Wannier eigenstates; an example of their structure 
(for v = 3) is shown in Fig. (1). For n = 1/3 one obtains 
a unique, fully saturated ferromagnetic ground state. For 
n < 1/3 only the Wannier states with a spatial over- 
lap have the same spin. The (degenerate) ground states 
hence consist of ferromagnetic clusters of arbitrary spin 
orientation. When the density is increased to n = 1/3 
the clusters touch and the degeneracy is lifted. Since 
Hjj\^g(N)) = and the kinetic part of H is diagonal 
in real space, the ground state © is localized. This is 
an explicit realization of Mielke-Tasaki's flat-band ferro- 
magnetism [l(| ■ Since the lowest single-electron band is 
flat only for 5 = tt/2 and is dispersive for 5 — 0, and 
the system therefore most probably conducting, we here 
encounter a metal-insulator transition as a function of 
magnetic field. 

Solution II. Correlated half-metal: The results of Ref. 
Q suggest that itinerant states are easier to realize at 
5 =/= tt/2. To investigate this point we analyze a group 
of solutions for flux 6 — tt, density n > 4/3, hopping 
tx, t\\ > and local potential e = l/(2tj_) — tx- Then 
the bare band structure consists of two upper flat bands 
Ei,E 2 , and a lower dispersive band ^(k) (Fig. 2b). 
For simplicity we discuss here only the case txtn =1/4, 
which implies E\ — E 2 = e + tx and £3 = 6 — tx — (1 — 
cos ak)/(2tx). Electrons in the dispersive single-electron 
band are created by the fermionic operators 



C2,k,, 



Cl.k.c 



k.cr 



\/2i?k 



yl — cos ak 
2£±i?k 



' c 3,k,, 



(7) 



with Rk = y/l + (1 — cos ak) 
canonical fermionic operators \li 



Defining non- 



localized ground states for densities n < 1/3, tx = tu =0 ^4i iCr = y/t\\ [ci )CT — Ci +BL}tJ — 2txe* 2 (ci+ ri)(T — Ci+ r2 , CT )] (8) 
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Eq. ([!]) is transformed into positive semidefinite form 



E 



ii. 



(9) 



(e + U + t x )N - N c [3U + At ±_ + 1 /t±] is the ground- 
state energy. For the positive semidefinite operator P = 

Sj=i (^id^J,! — Ai-T — Ai-l + 1) m ^1- ® to assume its 
minimum eigenvalue there must be at least one electron 
at each site. For N — 4N C electrons (n = 4/3) the ground 
state then has the form 

|*£'(4JV C )) = c[Y[A{_Jl}FnO) (lO.a) 

i 

= n^t-n<k,jio), (io.b) 

k s=l 

where F% creates two electrons with spin a on arbi- 
trary sites of each unit cell; c is a normalization factor. 
Since A[ creates one more spin a electron in each unit 
cell, every lattice site is occupied by a spin a electron. 
The electrons with spin — a are spatially extended, but 
they are localized in the thermodynamic limit. This is 
inferred from the long-distance behavior of the ground- 
state expectation value of the hopping term T r ^ a = 



-H.c). With j = i+rx, r/a = m and N c 



oo one finds r m ,_ CT = [(-l) m exp(-m/£_ CT )]/ y/1 + TJtZ. 
The one-particle localization length 



= -{ln[l + (2t 1 _) 2 - (2t ± ) v /(2^)2 + 2]}- 1 . (11) 

is finite; it increases almost linearly with l/t±. For N = 
4N C + AN electrons the ground state is given by 



AN 

n 



it 



,k„,-c 



(12) 



where n a can take the values s = 1,2,3. Since the op- 



erators c. 



t 

n a ,k Q 



add plane wave-type states to 



the ground state now contains genuinely itinerant spin 
—a electrons 1J]. This, together with the fact that 
5fi(N) = for 1 < AN < N c implies that the ground 
state l^l 1 ) is conducting. The net magnetization de- 
creases linearly with increasing electron density, l^g 1 ) re- 
mains the ground state up to A7V = N c (n — 5/3), where 
it becomes non- magnetic. For densities 4/3 < n < 5/3 
the ground state therefore corresponds to a correlation 
induced half-metal. Namely, while the 3N C electrons with 
spin a are completely immobile and iV c electrons with 
spin —a are confined to their localized Wannicr function, 
only the AiV — a electrons are itinerant, leading to a low 
carrier-density metallic behavior with a low spin polariza- 
tion. Since the conduction through this correlated half- 
metal involves only electrons of one spin species, such a 
system may, in principle, serve as a spin valve device, if 
contacted by metallic reservoirs 2ilZ|- Bearing in mind 



that at U — the electronic states at the Fermi level are 
dispersionless and hence localized, we find that a finite, 
repulsive Hubbard interaction can induce a localization- 
delocalization transition towards a half-metal. This is an 
explicit example for the correlation-induced conducting 
state conjectured by Vidal et al. Q. At fixed magnetic 
field it can be realized by tuning the sublattice poten- 
tial e. Similar solutions at S — tt can be deduced also 
at ii|i_L < 1/4, when the two upper flat bands are non- 
degenerate. 

Solution III. We now construct exact ground states 
for more general values of the magnetic flux, 5 G 
(— 7r/2, +7r/2), including zero flux, for electron densities 
n > 5/3. In particular, we will show that, by switching on 
a magnetic field, a non-magnetic ground state may turn 
into a non-saturated ferromagnet. We note that this oc- 
curs in the absence of any Zeeman coupling to the spin 
and is only due to the Peierls factor in the kinetic energy. 
Specifically we select the parameters t\\ = 0, t± < 0, 
6 = — cosS/t±, e = 6 — 6 _1 . In the non-interacting 
case one obtains a band structure with one upper flat 
band at E\ = e + 1 /6, and two lower dispersive bands at 



E- 



2.3 



e - 6/2 ± (b 2 /4 + t 2 ± +cos<5cosafc) 1/2 (Fig. 2c). 



Defining the non-canonical fermionic operators [li 



1 I 6c i+ a i0 . - c i+I . 2ja -e l 2 - c 1+ri!0 -e ^ 
\/2b j 6c i:CT - c i+r2j(T e~ l 5 - c i+riia -e l 5 



the ground-state energy E, 



in 



(13) 

(U+b)N-N c (3U+A/b+ 
26) and the Hamiltonian Eq. |T]) in positive semidefinite 
form are obtained as 



H 



^2 A*,i,gAL.i., 

a— ± \,<j 



UP + E, 



in 



(14) 



III. a Field induced localized jerromagnetism at n = 
5/3: For N = 5N C the ground state has the general form 



\* I J I (5N C ,5))=&E\5)\0) 



(15) 



where # = G\g\, with & a = U a =±Ui^ a ,i^ inserts 
two electrons with spin a and two electrons with spin 
— a in each unit cell. The first term of the Hamiltonian 
(fl"4|) therefore annihilates the ground state (fT5|) . In or- 
der to fulfill also PG^E^(S)\0) = 0, the operator &(5) 
requires a special form, which depends on the proper- 
ties of &: For zero flux (5 = 0) & creates a double 
occupancy on each site of the s — 3 sublattice, where 
the diamonds touch, thereby blocking any electron mo- 
tion along the chain. The operator E^(S = 0), which 
adds one more electron to each unit cell, can there- 
fore place electrons only on sublattices s = 1,2, i.e., 

^ f (°) = ni( a i4+r 1)ffM +A £ i+r 3 _, - i , 2 ) with arbitrary nu- 
merical coefficients ai,0; and spins oi l, a- h 2- Altogether 
the ground state is localized and non-magnetic with a 
high spin degeneracy and has the form 

[ni(t+^#(o)io)-(i6) 



\* I g II (5N c ,0)) = 



4 



For finite flux (8 ^ 0) & no longer creates a double occu- 
pancy on each site of the s = 3 sublattice, thus allowing 
electrons to move. To provide eigenvalue zero for the first 
two terms of H, (fT"4")) . E' (8 ^ 0) must be chosen such that 
it introduces one electron with fixed spin (say, j) on each 
site of the s = 3 sublattice. As a consequence there will 
be one t spin electron on all sites, i.e., these electrons are 
immobile (iV T = 3N c ,N l = 2N C ). For 5 ^ the ground 
state for N = 5N C is therefore a non-saturated ferromag- 
net with magnetization M oc (1 — AN/N C ) independent 
of S and a finite one-particle localization length £j (the 
proof proceeds as in [14j|). 

IILb Field induced itinerant ferromagnetism and 
insulator-metal transition at 5/3 < n < 2: The prop- 
erties of the ground state at S = remain valid even 
for N = 5N C + AN with 1 < AiV < N c since the 
ground state is now \# g (N, 6 = 0)) = [11^(^4+^,^ + 
A'4+ r ,,<T i2 )l*s( 5iV c,0)). The product is taken over AN 
arbitrary sites; the coefficients a[, (3[ and spins c^i, o\ i 
are arbitrary. Consequently, at 6 — one obtains again 
a localized ground state. For finite magnetic fields the 
ground state is instead 



\9 I "(N>5N C ,S)) 



AN 

n 



,k Q ,-( 



\9l"(5N e ,6)).{lT) 



Since the ground state now contains plane wave-type 
contributions and 5fi(5N c + AN) = for AN > 1, 
\^g II (N > 5N C ,S)) corresponds to a conducting, non- 
saturated ferromagnetic state with only mobile spin —a 
electrons. Hence the magnetic field induces an insulator- 
metal transition. In the metallic state the net magnetic 
moment decreases linearly with increasing density. 

Solution III therefore has the following properties: At 
zero magnetic field, in = 0, and a sublattice poten- 



tial 



tj_ - t 



_L ' 



but otherwise arbitrary t± < 0, 



U > 0, it represents a localized non-magnetic ground 
state over a continuous range of densities n > 5/3. 
By contrast, at finite magnetic field and the potential 
e = t±/ cos S — cos S/t±, but otherwise for the same pa- 
rameters, a non-saturated, ferromagnetic ground state is 
obtained. This state is localized at n = 5/3, but gapless 
for n > 5/3. For the latter density the majority spin 
electrons are immobile and only the minority spin elec- 
trons are itinerant. Therefore, by varying the magnetic 
field and the sublattice potential one can tune from a lo- 
calized, non-magnetic ground state in the density range 
n > 5/3 to a non-saturated ferromagnet, which is insu- 
lating at n = 5/3 and gapless for n > 5/3. 

In summary, by constructing exact ground states of 
electrons on the diamond Hubbard chain in a magnetic 
field we showed that this one-dimensional structure dis- 
plays remarkably complex physical properties which orig- 



inate from flat single-electron bands. The selected solu- 
tions describe flat-band ferromagnetism, correlated half- 
metal behavior with spin-valve features, and insulator- 
metal transitions. These properties do not depend on the 
Zeeman interaction [ls| . The virtue of tuning fundamen- 
tally different ground states through external magnetic 
fields or site-selective potentials thereby points to new 
possibilities for the design of electronic devices, which 
can switch between insulating or conducting and ferro- 
magnetic or non-magnetic states. 
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